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Critical exponents and unusual properties of the broken phase 
in the 3d-RP^ antiferromagnetic model.* 

H.G. Ballesteros, L.A. Fernandez, V. Martm-Mayor and A. Munoz Sudupe 

Depto. de Fisica Teorica I, Fac. de CC. Fisicas, Univ. Complutense de Madrid, 28040 Madrid, Spain. 

We present the results of a Monte Carlo simulation of the antiferromagnetic RP^ model in three dimensions. 
With finite-size scaling techniques we accurately measure the critical exponents and compare them with those of 
0{N) models. We are able to parameterize the corrections-to-scaling. The symmetry properties of the broken 
phase are also studied. 
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1. Introduction. 

AntiFerromagnctic (AF) models exhibit inter- 
esting characteristic properties. For instance, AF 
couplings on non-bipartite lattices, as the trian- 
gular Ising model, produce frustration. On bipar- 
tite lattices they can produce disordered Ground 
States, as for the AF three states Potts Model in 
three dimensions, which belongs to the XY Uni- 
versality Class. It is then possible to find new 
Spontaneous Symmetry Breaking (SSB) patterns, 
and maybe new Universality Classes. Our inter- 
est in these models is twofold: they apply to Con- 
densed Matter systems (spin glasses, helical and 
canted spin systems, ^He superfluid phase transi- 
tion, etc.) and four dimensional Field Theory, 
as they might offer some insight on the formula- 
tion of non asymptotically-free interacting theo- 
ries. For the AF RP^ model, the finite -size scaling 
analysis close to the transition, strongly suggests 
that the action's 0(3)-symmetry is broken, yield- 
ing an S0(3)/{1} SSB pattern. The pcrturbative 
prediction for such an SSB pattern is that the 
transition can either be on the 0(4) Universality 
Class or that it should be first order or tricriti- 
cal0. 

2. The model. 

We consider a system of three components, nor- 
malized vectors, placed on the nodes of a cubic 



lattice. They interact through a first neighbors 
coupling 
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This action presents a global 0(3) symmetry, 
and a local Z2 one {cTi —<Ti). From Elitzur's 
theorem, it follows that we are really studying 
RP^ variables. For positive /?, this model suffers 
a first order phase transition, all spins aligned or 
anti-aligned with an arbitrary direction. This cor- 
responds to the nematic phase of liquid crystals. 
At /3 w —2.41 the system undergoes a second or- 
der phase transition |^,^ . The symmetry descrip- 
tion is however, more complicated. Let us call 
the spin placed on {x, y, z) even or odd, accord- 
ing with the parity of x -\-y + z. At /3 = —00, the 
ground state consists on all the spins on the, for 
example, even sublattice, aligned or anti-aligned 
with an arbitrary direction, while odd spins lie 
randomly on the perpendicular plane. This state 
has a remaining 0(2)-symmetry, but the sym- 
metry between even-odd sublattices is broken. 
However, thermal fluctuations induce an interac- 
tion between the spins on the plane sublattice, as 
an alignment allows stronger fluctuations on the 
other sublattice. 

As the natural variable is the tensorial product 
of (T by itself, we consider the traceless tensorial 
field 



Tf^ = apdf 



(2) 
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We define the normalized non-staggered and 
staggered magnetizations as 



M 



V ^ 



{x,y,z)y 



x,y,z 
x,y,z 



In the simulation we actually measure 



M = (^VtrM2^ ^ 

and the associated susceptibility 
X = V(trM2) . 



(3) 
(4) 

(5) 
(6) 



We have found very useful the measure of the 
second momentum correlation length defined as 



X/F-I 
4sin2(7r/L) 



1/2 



(7) 



where F is the Fourier transform of the correla- 
tion function at minimal non zero momentum, as 
well as the scaling function (related to Binder's 
parameter) 



((trM^ 



(trM2) 
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The staggered counterparts of the quantities 
(^H) are defined analogously. 

We have also measured the first and second 
neighbors energies: 



^1 ~ ~ E ('^» ■ '^i)^' 



3V ^ 
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i?i has been used for the Ferrenberg-Swendsen 
extrapolation method, and for calculating /3- 
derivatives. i?2 is useful to obtain information 
about the structure of the broken phase. 

3. The simulation. 

We have used a Metropolis update. Cluster 
methods are feasible but not efficient. 

We have measure the integrated and exponen- 
tial autocorrelation times for x ^^id Xs,- The 



integrated times satisfy r^'* < r™* < r™* (for 
instance, in the larger lattice the ratio is 1:3:7). 
However, the exponential times are almost equal, 
and they are near t™*. So, we have confidence 
that the latter is the larger autocorrelation time. 
Further details can be found in refs. 

In table || we present the number of Monte 
Carlo iterations performed at each lattice size, 
together with the corresponding autocorrelation 
time. The expected behavior r c>c jg -v^ell sat- 
isfied. 



Table 1 

Number of Monte Carlo sweeps performed for dif- 
ferent lattice sizes. Measures have been taken ev- 
ery 10 sweeps. We have discarded in each case 
about 200t'"' iterations for thermalization. 



L 


MC sweeps(xlO'') 


^int 


# of 


6 


6.71 


7.37(3) 


910,000 


8 


17.07 


11.41(4) 


1,496,000 


12 


6.51 


24.9(2) 


261,000 


16 


22.14 


44.5(3) 


498,000 


24 


8.77 


107(3) 


82,000 


32 


28.51 


175(6) 


163,000 


48 


3.93 


410(20) 


9,600 



4. Finite size scaling techniques. 

For an operator, O, with critical exponent xq-, 
we have 

(0(L, /?)) = L^-l'^ {Fom, P)/L) + 0{L--)) . 

Order ^(oo, j3) ^ terms, are negligible in the crit- 
ical region. 

Measuring at two lattice sizes L and sL and 
computing the quotient 



Qo 



{0{-^L,(i)) 



(11) 



we can eliminate the scaling function Fq just by 
choosing the f3 value such that the ratio between 
correlation lengths is s, obtaining 



(12) 
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5. Exponents 

To compute the v exponent we consider op- 
erators like the j3 derivatives of the correlation 
length or the logarithm of the magnetization 
{^Ai/dP = v + l, XdiogAZ/d/s = !)■ In table | 
we present in the second column the results ob- 
tained from Qdjs/d/sl " 2. The values using 
other operators involving ^ are similar but slightly 
worse. In the case of operators involving magneti- 
zations, the statistical errors are smaller, but the 
finite-size effects are greater. We do not find any 
significant deviation between the staggered and 
non-staggered channels what supports the equal- 
ity V = Vs necessary to define a single continuum 
limit. 

In the case of magnetic exponents, we expect 
a different behavior for the staggered and non- 
staggered channels. We obtain the respective 
77s and 77 exponents from 7(s)/i^ or I3(s)/v us- 
ing the scaling relations 77(3) = 2 — 7(s)/^' or 
77(s) = 2 — 0? -t- 2/3(s) / V. The quality of the results 
depends on the observable measured. With a 
good selection of both the O operator and the 
definition of correlation length the corrections-to- 
scaling terms can be largely reduced. We should 
remark that the operator used to measure the cor- 
relation length in practice can be any quantity 
that scales linearly with L at the critical point, 
such as kL. However, the more interesting prop- 
erty of the method we use is that even if the op- 
erator O is a rapidly varying function of the cou- 
pling at the critical point, as the magnetization 
or the susceptibility are, the statistical correlation 
between Qo and allows a very precise measure 
of the critical exponents. 

In table ||, (col. 3), we display the results for r]s 
using Xs as operator and as correlation length. 
Regarding 77 exponent the use of ^ or as corre- 
lation length produce large corrections-to-scaling. 
We have realized that those effects are strongly 
reduced when using kL as the correlation length 
operator (see col. 4). 



6. Scaling corrections 

At very large volume, scaling functions for dif- 
ferent lattices, such as S^s/L, S,/L, k or Kg should 



Table 2 

h', rjs and 77 exponents obtained from ([l^ ) using 
as operators d^s/d/3, Xs and x respectively. As 
correlation length we take ^s in all cases but the 
latter where we have used kL. The pairs are of 
type (L,2L). 



L 


J/ 




V 


6 


0.786(6) 


0.0431(10) 


1.332(6) 


8 


0.785(4) 


0.0375(7) 


1.332(4) 


12 


0.789(8) 


0.0357(17) 


1.321(13) 


16 


0.786(6) 


0.0375(12) 


1.340(5) 


24 


0.77(2) 


0.038(5) 


1.334(18) 



cross at /3'^. Scaling corrections shift the cross- 
ing point, corresponding to a pair of lattices 
(Li,L2 = sLi), an amount given by 



Li ,L2 



1 



1 



■L7 



(13) 



As the crossing point of correlation lengths and 
Binder-like parameters tend to the critical point 
from opposite sides, it is sensible to perform a 
global fit, with the full covariance matrix to take 
into account the statistical correlation of all the 
data. We have fitted our data either fixing the 
small lattice, Li, or an s value. For instance, for 
Li = 8, we obtain 



xVd-o.f. 



-2.4085(3) , 
0.86(4) , 
12.2/14 . 



(14) 



Using also the Li = 6 data or fixing s — 2 the 
results are compatible (see ref. [^). 

After measuring the, universal, corrections-to- 
scaling exponent uj we can estimate the finite-size 
effects on the critical exponents using 



X 
V 



(2L,L) 



cx L 



(15) 



From the results quoted in table y we conclude 
that the values for v are rather stable when in- 
creasing the lattice size, but a proper accounting 
of finite-size effects is crucial in the case of 77 ex- 
ponents for some operators (see also ref. Q where 
this method is applied to three dimensional O(A^) 
models). 
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Anyway, the presence of an L dependence, even 
when it is not measurable, requires an increas- 
ing of the error bars for a safe determination of 
the systematic error due to finite-size effects. We 
summarize the results giving the value for the 
(16,32) pair, with a second error bar that cor- 
responds to the increasing of the error due to the 
infinite volume extrapolation. 



V = 0.783(5)(6) , 

= 0.0380(12)(14) 
•q = 1.339(5)(5) . 



(16) 



Some minor differences between the values in ta- 
ble H and (^ exist because here we average be- 
tween the results of several operators not dis- 
played in the table. 

The exponent u is two standard deviations 
apart from the value for the three-dimensional 
0(4) model {vo(i) = 0.755(8)[D) what supports 
that the AF transition in the RP^ model belongs 
to a new Universality Class, but the opposite can- 
not be completely discarded. 

7. Vacuum structure 

In the P = —oo limit, there is a breakdown of 
the symmetry between the odd and even sublat- 
tices, but a global 0(2) symmetry still remains. 
At the critical region, the presence of fluctuations 
can change this picture. 

To study the structure of the broken phase near 
the critical point we have analyzed the FSS of two 
operators 

A = ((trM,M)2), 
B = (tr[M3,M][[M„M]t). 



(17) 
(18) 



The A operator depends on the difference of 
the magnetizations squared of both sublattices 



(trM,^„ 

try were not broken it should scale as L 
obtained value 



trM^jj^), so, if the even-odd symme- 

The 



xa/i^ = -3.389(15) 



(19) 



agrees with — 2(/3 + fis)/v and discards a remain- 
ing even-odd symmetry. 

The B operator measures if the tensors M and 
Ms commute (what is equivalent to the commu- 
tation of Mcvcn with Modd)- If they did not, a 



simultaneous diagonalization is not possible, and 
a remaining 0(2) symmetry is discarded. The ex- 
pected FSS for an unbroken 0(2) symmetry cor- 
respond to xb/v = -d - 2j3s/v = 4.04(2) while 
the value obtained numerically is 



xb/v = -3.406(11) 



(20) 



To check if both breakdowns occur deep in the 
broken phase, we have studied the second neigh- 
bor energy E2 for each sublattice, as well as the 
eigenvalues of the magnetization tensors. 

The analysis of the E2 histograms show an un- 
ambiguous signature of breakdown of the even- 
odd symmetry along all the broken phase. It is 
interesting to remark that the energy of the less 
aligned sublattice takes for fi — > —00 the value 
0.532(5) to be compared with the (3 — —00 value, 
0.5. The presence of fluctuations acts as an effec- 
tive coupling that favors an alignment. 

The intensity of that induced interaction in- 
creases when approaching to the critical point. 
To know where it actually produces a breakdown 
of the 0(2) symmetry, we have studied the eigen- 
values of the M matrix. We have observed that 
the difference between the two smaller eigenval- 
ues goes to zero when the lattice size grows, for 
/3<— 4. At (3 — —Z a. nonzero value is not com- 
pletely excluded but we would need lattice sizes 
larger than 48 to obtain a reliable measure. In 
order to clarify this point, we project to study an 
extended model where an explicit second neigh- 
bor coupling is added. 
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